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^ ' Abstract 

O , In terms of proposed by authors general-relativistic kinetic model of 

■ baryon production in expanding primordially symmetrical hot Universe 

1 ' calculates distribution function of extra-massive bosons and concerned 
, with it variables. 

CJ ! 1 Introduction 

cr. 

i In previous paper [I] in terms of theory , developed earlier by one of the authors 
was proposed an accurate kinetic model of baryon production in primordially 

I baryon-symmetrical Universe on basis of field model of interaction of particles of 

T-H ■ type SU(5). Here and further references to formulas of authors previous paper 

' [T] are indicate as (I,N) , where N - number of formula of stated paper. 

m ; 

(N . 

^ ; 2 Weak violence of charge symmetry in stan- 
: dard SU(5) model 

I Using weakness of violence of CP-symmetry 

Tl ; Ar = r - f < 10~^ 

> ■ 

. and thereby smallness of chemical potentials (Aa ^ Ar ^ 1), we can reduce 

' referred above system of intcgro-difFcrential equations to a system of linear 

. integro-differential equations. In present paper we will consider simple model of 

" " ' baryogenesis, in which CP-invariance is violated in decay channels of only one 

type of bosons, A-boson. In this case equations (1.40), (1.42) assume following 

form: 
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Here separately wrote out equations for anoquarks- ([2]), and katoquarks 
At derivation of these equation was allowed that: 

Pix,x)^2p{x). (4) 

Let us write out formulas for coefficients that belong to these equations in 
linear in Ar, A approximation. In zero approximation: 

m^MP,t), (5) 

where: 

« m - ^ In 1 + exp(p+/r) 

^) - p exp(p/2r)[l + exp(p_/T)] 

and it is necessary to substitute: 
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Further: 



where: 
and 



/^(^) = /^(0;P,^)+^ exp(gj,/r) 

^xv^7 , , y [exp(£;x/r) - 1]' ^ ^ 

/°(0;P,i) = [exp(i?x/T)-l]-i (9) 
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Integrating kinetic equations for X-bosons ([1]) with initial conditions (1.43) 
in linear in Ar, A approximation we will determine: 

/x(P, t) = /^(O; P, 0) + <5/o(P, t) + <5/(A), (10) 
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where Sf{X) - linear in A functional, hence: 

6{-X) ^ -S{X), (11) 

and it doesn't contribute in final output; and 6f - deviation from equilibrium 
of distribution function: 

t 

Sf{P,t) = -e-*(^'*) / e*(^'*''/o(0;P,i')^^i' (12) 





- deviation from equilibrium of boson distribution function in symmetrical plasma 
A = 0, and incorporated notation: 

X 



To J ^a2(t') +p2/m 







To - proper time of free X-boson's decay: 

To = — -^{mxa) \ (14) 

Since fermions lay in condition of and are ultrarelativistic, their concentra- 
tions are equal: 

oo 

Af ^ ^ [ , , ^ P^dP. (15) 

7r2 J exp(-A + P/r) - 1 ^ ' 



In conditions of weak violence of CP - invariance, when Aa <C 1 we will receive 
from here approximately, separating at degrees of smallness A : 

3r3c(3) r3 

2 TT^ 6 

ANa « Aa^. (16) 

In standard SU(5) model probabilities of X-boson's decay in lepton, {qe), 
and quark, (qq), channels are the same[l, i.e.: 

l-2r = 0i(Ar). (17) 

Subject to this factor and relation (1161) in standard SU(5) model from equations 
for anoquarks- ©, and katoquarks -([3]) we can receive one closed first-order 
equation on variable B = u + 2d, i.e., on excessive concentration of baryons: 
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^subject to colors and charms 
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= ^/p^W/orfP, (18) 



integrating which subject to initial conditions (1.43), we will obtain: 

oo / oo \ 

AAAb(oo) = y" exp I - y" ^{t')dt' j G{t)dt, (19) 



where 



oo 

*(0 = ^ 1 P'$/o/3orfP, (20) 





oo 

1 ' T)2 



G(t) = ;^ y P^$(5/dP. (21) 



Now it remains to calculate variable Ss = Sng/S, where entropy density of 
ultrarelativistic gas is equal: 

S = ^NT^ ^ 5 - —NT\ (22) 
45 45 

By that the task is formally solved. We should note that in contrast to 
papers in which distribution function of X-bosons was modelled by quasi- 

hydrodynamic distribution and was obtained by method of numerical integration 
of kinetic equations, distribution function of X-bosons in present article obtains 
by strict integration of kinetic equations and B(po) determines in quadratures: 

- (23) 

3 Analysis of solution 

Further we will assume that firstly X-bosons decay in general on intermediate 
stages of expansion, when T ^ mx, and, secondly that part of X-bosons suf- 
ficiently small in comparison with general number of particles: Nx <C N, so 
with enough accuracy grade we can assume cosmological plasma on decay stage 
of X-bosons is ultrarelativistic. This gives us laws of variation of dimensioned 
factor and temperature in due course: 

a{t)^aoVi; T = ro^^, (24) 



4 



where To (• (1-44)) is plasma's temperature in Planck miits on Planck moment 
of time: 

■ 45 



Let us explore obtained in previous section linear solution for that we will 
come over from variables t,p to new variables ry,^: 

t - Tom (26) 



—pmx^; 



mx£,y/n) 



(27) 



where tq - proper time of free X-boson's decay (jl4p . in that way to valuation 
of dimensionless time 77 = 1 it corresponds time t = tq. Let us incorporate 
dimensionless parameter a , depending from constants of field theory q 



T{ro) 



To 



(28) 



which is equal to relation of X-boson's mass to temperature at a point of it's half- 
decay, where x - dimensionless parameter of 1 order, depending from parameters 
of field theory, a - coupling constant ( 0.1 0.01). 
Then we will receive expressions: 
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/o(t,P)= l + eM^^i^ + e) 



/3o(P,i) = /3o(e,?7,'T) 



= — In 



l + cxp((T/2(y^TF + ^)) 



l + exp(a/2(v/;^TF-0) expK/2) 



(29) 
(30) 

(31) 



On fig. [T]is shown relation of statistical factor /3o(^, ?7, <t) pip from parameters 
cr ^ =_p. 



^Incorporated in ^ parameter a is equal to quadrate of cr, which uses in present article. 



5 




3.1 Boltzmann distribution 

Let us calculate function ^{F,t) p^ . going over variables £,,7] according to 
formulas CT. (P7|) : 

mt) = mv)=J^^^^^^^. (32) 

In Boltzmann approximation, substituting statistical factor /?□ ^ and executing 
integration we will receive: 
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3.2 Deviation of X-bosons from equilibrium 

Substituting function $(^, i]) in form ((33)) in expression p2|) for X-boson's distri- 
bution function deviation from equilibrium and going over variables p6p . ()27p . 
we will obtain: 

1 7 exp($(e,r,')+av/7+e^)rf^' 



6 



On Fig. [21 [3] results of numerical integration of expression ([M)) for relative 
deviation of distribution function from equilibrium, Sf/fO, subject to time, 
77, and momentum, ^ for value a — 1 are shown. Let us notice that fermi- 
distribution is substituted to boltzmann distribution only in expression (j32[) for 
function $(^, 77)). 




Fig. 2. Relative deviation of distribution function of X-bosons from equilib- 
rium (5/(^, r], cr)//o under a — I. 

From dduced results follows that: 

1. Maximum of relative deviation of distribution function from equilibrium , 
((5///o), falls at lesser values of momentum coordinate ^ < 1; 

2. On initial stages relative deviation of distribution function from equilib- 
rium grows with time 77; 

3. However at achievement certain critical sufficiently great value of time 
relative deviation of distribution function from equilibrium begins to de- 
crease, - it associates with a fact, that by this moment of time X-bosons 
start to decay; 

4. At very great times weak maximum of relative deviation of distribution 
function from equilibrium appears and starts to migrate to the range of 
great values of momentum coordinate. 
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Fig. 3. Evolution of relative deviation of distribution function of X-bosons from 
equilibrium 

Sf{^, rj, a) I fo under a — 1. Thin line - rj — 0.3, normal line - rj ~ 1, heavy line 
- ri — 3, thin dotted line - rj = 10, heavy dotted line - r/ = 20. 

Let us notice that specific features of behavior of distribution function's 
deviation were exhibit but were not explored in paper ^2^ . It is necessary to allow 
that under great values of time coordinate r/, absolute magnitude of density of 
X-bosons number becomes vanishingly small. On Fig. |3] - HI are shown results 
of numeric integration of expression (j34p for density of X-bosons deviation from 
equilibrium, £,^Sf, subject to time, r], and momentum, ^ for values a ~ 0.3, 1, 3. 
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Fig. 4. Density of non-equilibrium X-bosons number dN = ^'^6f{^, r], a) under 
£7 = 0,3. 
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Fig. 6. Density of non-equilibrium X-bosons number dN = £,^5f{£,, r], a) under 
cr = 3. 

On Fig. [T]- 19] are shown results of numeric integration for evolution of X- 
bosons density's deviation from equilibrium, £,^Sf, calculated by formula (p4)) . 
for values a = 3,1, 0.3. 
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under a = 0.3. Thin line -r] = 0.3, medium line - rj = 1, heavy line - r] = 3, the 
heaviest line - r] = 10, dotted line - t] = 30. 
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Fig. 8. Evolution of number density of non-equilibrium X-bosons dN = S,^6f{^, r], a) 
under a = 1. Thin line -r] = 0.3, medium line - r] = 1, heavy line - r] = 3, the 
heaviest line - r] = 10. 
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Fig. 9. Evolution of number density of non-equilibrium X-bosons dN = £,^Sf{£,, rj 
under = 3. Thin line -rj — 0.3, medium line - rj = 1, lieavy line - r/ = 3. 



3.3 Exponential factor in formula (llOh 



Further from formulas (rr9l) - (|22l) it is seen that final solution defines not by 
function itself $(P, t), but by it's time derivative, calculating which we will find, 
proceeding to variables ^, ry according to formulas (f26| . ([ST)) : 



i.(P,0 = -^^^ai2l^ (35) 
or, in Boltzmann approximation: 

Thus, calculating an integral by momentums we will receive expression for func- 
tion ^1/(77, ct) (|20p in Boltzmann approximation: 

vl/(,7,a) = ^^^i^i(ay^), (37) 

where Kn{z) - modified Bessel function (Bessel imagunary function) [?]: 

00 

K^{z) ^ J e-'''°^^''coshvudu, (38) 
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at that following recurrence equations are correct: 



(39) 



Thus, exponential factor in formula (|T9)) is equal: 

oo / 

"^^{1], a) = exp(- J a)dt') = exp 



2Nx 



V 



Considering that: 



37r 



(40) 



(41) 



it is not hard to receive limitary relation: 



^'e(?7,(T) K. exp 



(^07-^0), 



(42) 



as well as asymptotic expansion of function 5'e(f7, c) at small values of argument 

CTy^: 

3iVx 2Nx(Trf/'^^ 



^el??, cr) ~ exp 



37r2 



(aV^^O). 



(43) 



In other limitary case accounting asymptotic expansion of Bessel functions 
at great values of argument we will receive: 

*e(?7,fT) w exp (^^^V2^7f/^e-''^^ 



1 - ^\/2^?7^/4e""v^ [a^ ^ c5o). 



(44) 



On Fig llOl dependence of exponential factor from variables 77, cr, received by 
numerical integration for Nx = 1, is shown. 



13 




oo 

Fig.lO. Function ^e{v, = exp(- / "i^dt') at Nx = 1. 

t 

On Fig. [Tl]is shown an evolution of exponential factor (77,(7) at various 
values of parameter ct, 
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Fig.ll. Function ^'e(77,cr) = exp(- / -^dt') at Nx = 1. Thin line: a = 0.3, 
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medium line: a = 1, dotted line: cr = 3. 

On Fig. [T2] influence of number of X-boson types on exponential factor, Nx is 
shown. 
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Fig. 12. Influence of number of X-boson types, Nx, on exponential factor 

oo 

\I/i(r;,cr) = exp(- / '^dt') at a ^ I. Thin line: Nx = 1, heavy line: Nx = 3, 
t 

dotted line: Nx = 12. 

On Fig. [T2] are shown vahies of exponential factor, calculated by numerical 
methods and using asymptotic estimates P5)) . 
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Fig. 13. Comparison of exponential factor values, calculated hy numerical metli- 
ods (thin line) and by asymptotic estimates if43)) -heavy line and - dotted 
line. 
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